Ground-State Fidelity and Bipartite Entanglement in the Bose-Hubbard model 
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We analyze the quantum phase transition in the Bose-Hubbard model borrowing two tools from 
quantum-information theory, i.e. the ground-state fidelity and entanglement measures. We consider 
systems at unitary filling comprising up to 50 sites and show for the first time that a finite-size scaling 
analysis of these quantities provides excellent estimates for the quantum critical point. We conclude that 
fidelity is particularly suited for revealing a quantum phase transition and pinning down the critical point 
thereof, while the success of entanglement measures depends on the mechanisms governing the transition. 



A few years ago some key works [J] initiated a new vein 
of research using concepts borrowed from Quantum Infor- 
mation Theory in the analysis of quantum phase transitions 
(QPT, i.e. phase transition driven by quantum as opposed 
to thermal fluctuations). The best known examples are 
no-doubt the measures of entanglement, which quantifies 
the strength of quantum correlations between subsystems 
of a compound system and represents a basic quantum- 
computational resource J2]. A more recent proposal is 
based on the fidelity, a key parameter in the characteriza- 
tion of the performance of logical quantum gates JD]. The 
main advantage of this tool lies in the fact that, being a 
purely Hilbert-space geometrical quantity, it does not re- 
quire any a priori knowledge of the correlations driving 
the QPT, or of the order parameter thereof 0]. 

While most of the works in this relatively new field fo- 
cus either on fermionic models JH 0, 0, Hi |9|] or on spin 
models Gl S ES O, E B that can be often effectively 
posed as free spinless fermionic systems, bosonic models 
went somewhat unaddressed so far. Two exceptions in this 
respect are Refs. ILHl and [151, which propose the study of 
the hallmark QPT of the Bose-Hubbard (BH) model using 
respectively entanglement and Loschmidt Echo, a quan- 
tity kindred to fidelity. We note that the latter also pro- 
vides an experimental scheme to measure fidelity. Also, 
the crossovers characteri zing the ground-state properties of 
the attractive BH model lfl6ll are investigated in terms of fi- 
delity in Ref. I117ll . This substantial lack of attention does 
not make justice of this paradigmatic bosonic model. In- 
deed, on the one hand the BH model has a clear experi- 
mental relevance, being standardly realized in terms of op- 
tically trapped ultracold atoms ltl8ll . On the other hand, it is 
a genuinely many-body model which in general cannot be 
reduced to an effective noninteracting model, hence pos- 
ing a significant computational challenge. These features 
make the BH model ideal grounds for investigating the ef- 
fectiveness of Quantum Information tools in the study of 
QPT. 

This is the aim of the present work. We focus two spe- 
cific topologies, i.e. the one dimensional lattice with pe- 
riodic boundary conditions (ring) and the completely con- 



nected graph (CCG), i.e. a model with the same hopping 
amplitude across any two sites. A twofold reason makes 
the latter a convenient benchmark [fsL 151 [Toi [Til [H 
First, the critical lines of its zero-temperature phase dia- 
gram are known analytically in the thermodynamic limit 
]2fj]; second, its high degree of symmetry allows for a sig- 
nificant reduction of the relevant Hilbert space. As to the 
ring, a very interesting recent proposal J21 ] turned it from 
a convenient theoretical idealization to an experimentally 
realistic system Jl6, 1_7, 22]. More specifically, we investi- 
gate the relation between the behaviour of the ground-state 
fidelity and bipartite mode entanglement — both in the di- 
rect and reciprocal space — and the superfluid-insulator 
QPT taking place in the pure BH model at inte ger fill- 
ing. An efficient use of the system symmetries Il23ll al- 
lows us to apply exact diagonalization algorithms to rings 
and CCG's comprising up to 12 and 50 sites, respectively, 
and containing an equal number of bosons. We also con- 
sider the Gutzwiller mean-field approximation to the BH 
model, whose phase-diagram is essentially the same as that 
of the CCG 120]] . Our conclusions can be summarized as 
follows. The finite-size scaling of the position of suitably 
chosen extrema of both the fidelity and the bipartite en- 
tanglement provides a satisfactory estimate of the critical 
point of the QPT, for both geometries. Consistent with 
what argued in the seminal Ref. |4|], fidelity turns out to 
be the candidate of choice when it comes to singling-out 
and quantitatively pinpointing the critical parameter of a 
possibly unknown QPT. Some of the appealing features 
that tip the scales in favour of fidelity in this task include: 
a very intuitive definition; a straightforward evaluation; a 
clear telltale of the transition; the robustness against mean- 
field approximation (note indeed that, unlike bipartite en- 
tanglement, the fidelity is defined also for the product trial 
state inherent in the Gutzwiller ansatz). As to the bipartite 
mode-entanglement, its more complex definition and less 
straightforward evaluation make it a less handy tool for lo- 
calizing the quantum critical points. Moreover, as it is often 
pointed out in the literature (see e.g. Ref. lfl2ll ) a system- 
independent recipe based on entanglement measures seems 
to be lacking. However, the interest of such measures also 
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relies in the insight that they can provide about the mech- 
anisms driving the (possibly otherwise identified) quantum 
phase transition [6]. For instance, the failure or success of 
the recipes considered in this work can be explained based 
on the features of the QPT characterizing the system under 
examination. 

The Hamiltonian of the pure Bose-Hubbard model (i.e. 
on-site interactions, no offset in the local potential) reads 
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where a m , a' m and n m = a) m a m are lattice boson oper- 
ators which destroy, create and count bosons at site m. 
As it was argued in Ref. ll24ll and later confirmed exper- 
imentally 1181 . Hamiltonian £[]) describes ultracold atoms 
trapped in an optical lattice and its only effective param- 
eter J > 0, i.e. the tunneling amplitude to boson (re- 
pulsive) interaction ratio, is directly related to tunable ex- 
perimental parameters such as the atomic scattering length 
and the intensity of the laser beams providing the optical 
confinement Il25ll . The adjacency matrix A describes the 
coordination of the M sites composing the lattice, being 
nonzero only for adjacent sites. As we mention, we con- 
sider two topologies, i.e. the ring and the CCG, whose 
adjacency matrices read respectively A rnm i = 8\ m _ m >\\, 
where \M — 1| = 1 owing to the periodic boundary con- 



ditions, and A Ti 



2(1 - 8 m , m >)/{M - 1). The nor- 



malization factor in the latter ensures that the generalized 



coordination number z m = ^m'=i A m m > equals 2 in both 
cases, so that in the thermodynamic limit of M — > oo the 
CCG is actually an infinite range mean-field approximation 
to the ID system ll20ll . Indeed the critical boundaries in the 
phase diagram of the CCG 11201 are the same as those of the 
mean-field approximation ensuing from the decoupling as- 
sumption a^a m / ~ a] n {a m ,) + a m >{a\^ - (a^)(a m /), 
where (•) denotes expectation on the ground state of the 
system 

JHStI]- For a homogeneous lattice like the ring, 
the resulting mean-field (or Gutzwiller) Hamiltonian is the 
sum of M identical on-site contributions which, dropping 
the site subscripts, reads 



u = m 



n(n — 1) — fin — Uaia + a*) 



(2) 



where a is the mean-field parameter that is related to 
the relevant ground-state 1 1 J / ) by the self-consistency con- 
straint a = (^/\a\*S>), and [i is the chemical potential, 
i.e. a Lagrange multiplier controlling the total number of 
bosons N = ^2 m n m . Note indeed that since [N, H] = 
^ [N,7i], the boson population is not a good quan- 
tum number in the mean field approximation. Conversely, 
the ground-state of Eq. CO can be studied within fixed- 
number Fock spaces. This reduces the in-principle in- 
finite Hilbert space of the system to the large but finite 
size V(M,N) = ( N+ ^~ 1 ). The phase diagram of the 



an extended compressible superfluid phase and a series 
of incompressible insulating Mott lobes. The latter corre- 
spond to commensurate populations, i.e. to integer fillings 
v = N/M, and their boundaries are given by the crit- 
ical lines \i = /x±(J) = ±[E±(J) — E C (J)], where E c 
and E± denote the ground-state energy relevant to the com- 
mensurate (N = I'M) and defect states (N = uM ± 1), 
respectively lEH l29ll . On ID system the transition at in- 
teger filling occurring at the tip of the lobe belongs to the 
particularly elusive Berezinskii-Kosterlitz-Thouless (BKT) 
class l20h . The precise location of the critical point thereof 
requires ad-hoc procedures that rely on an a priori knowl- 
edge of the mechanisms driving the transition lE^ll . One of 
the best estimates for v = 1, = 0.26 ± 0.01, has been 
obtained via numerical st rong -coupling expansions of jjl± 
of remarkably high order I28II . As we mention, the bound- 
aries of the mean-field Mott lobes are known analytically 
ll20ll . In particular, the tip of the v = 1 Mott lobe is at 
J x = 3/2 - y/2. 

The critical point can be also studied considering in- 
teger fillings only, v = N/M £ N. Indeed, order param- 
eters such as the superfluid or condensate fraction feature 
extrema in the suitable derivative with respect to J [f30h 
that are expected to diverge in the thermodynamic limit. A 
finite-size scaling of the location J M of such an extremum 
is expected to result in the critical point J^, 



M = f(J M ) = C 1 \J M -J a 



C u ri>0. (3) 



We checked this expectation and our results, which we do 
not show here, actually provide estimates of that agree 
very well with the best known estimate for the ID lattice 
1 280 and the analytic result for the CCG. 

As we mention, in this work we are interested in the in- 
sight provided by observables borrowed from Quantum In- 
formation Theory, i.e. the fidelity and the bipartite mode- 
entanglement. The former quantity has a remarkably sim- 
ple definition, being nothing but the modulus of the over- 
lap of two ground states relative to two different choices 
of the Hamiltonian parameters J T ( J, J') = |(\I/jj\I/jv)|. 
On finite size systems, a drop in the fidelity corresponding 
to two arbitrarily close parameter choices J' = J + 5 J 
is expected to signal a precursor of the QPT 0, Hi |9fl. A 
perhaps more effective indicator is provided by the peak 
in the "density" of the second derivative of the fidelity 
<S(J) = lim 5J ^ 2[1 - T(J, J + 5J)]/{M5J 2 ) §]. The 
location J M of such peak on a size-Af system signals the 
QPT provided that the S(Jm) grows with increasing size, 
thus resulting in a divergence in the M — > oo limit. 

The bipartite mode-entanglement (henceforth simply 
entanglement) is the von Neumann entropy 8 (J) = 
— Tr[p*( J) In p* (J)]/ In I?*, where p* is the reduced den- 
sity matrix corresponding to one of the subsystems inher- 
ent in a system bipartition, and D* is the size of the rel- 
evant Hilbert space. We consider both the spatial modes 



BH model, usually drawn in the \i — J plane, comprises (SM) a rn and the quasimomentum modes (QM) b q 



3 



Y2 m e l ^ mq a m . In both cases p* is obtained by trac- 
ing out of the density matrix p{J) = \^Sj){^fj \ the degrees 
of freedom of all but one mode, so that T>* = N+l. In the 
first case we choose one of the (equivalent) spatial bosonic 
modes. In the second case the untraced degrees of freedom 
correspond to the mode q = 0, that describes the system 
ground-state in the noninteracting limit. 

We now illustrate our results, which have been obtained 
by numerically diagonalizing Eq. (TJ. An efficient use of 
the Hamiltonian symmetries allows us to consider rings 
(CCGs) comprising up to M = 12 (M = 50) sites at 
unitary filling. We remark that "D(50, 50) m 10 29 . Fig. 
Q] shows the Mott lobe corresponding to filling v = 1 for 
both the ring and CCG. The filled regions represent the re- 
sults in the thermodynamic limit as reported in Refs. ]28|] 
and lEM E^l . respectively. The solid lines refer to finite- 
size results as described by the color code, which applies 
also to the subsequent figures. Fig. [2] shows the behaviour 
of the fidelity J 7 (J, J + 5 J) and of the relevant fidelity 
derivative density <S(J) JI i, ED] for the ring and CCG. 
The extrema of such quantities feature a scaling behaviour 
that is consistent with the hypothesis that they correspond 
to precursors of the quantum critical point. This is sig- 
nalled by a shaded stripe (corresponding to the estimate in 
Ref. [28]) and a vertical line for the ring and the CCG, 
respectively. The same applies in Figs. [3] and |4] Fitting 
the locations Jm of these extrema as in Eq. © results into 
Joe = 0.257 ± 0.001 = 0.086 ± 0.005) for the ring 
(CCG) which agrees very satisfactorily with the known re- 
sult J x = 0.26 ± 0.01 H (</oo ~ 0.0858). This is 
illustrated by the black plots in Fig. |U We also checked 
the extensivity of the peak of the fidelity derivative, mea- 
sured by the exponent 7 in the fit S{Jm) ~ C 2 M 7 . This 
is in general expected to be related with universal quan- 
tities and, for the previously considered (effectively free) 
models, it equals 1 0, 0, [32]. In the case of the ring we 
find a slightly but definitely superextensive (7 > 0) be- 
haviour, 7 = 0.087 ± 0.009. In the case of the CCG, 
7 = 0.749 ± 0.007 signals a markedly superextensive be- 
haviour. 

As it is clear from Fig. [3j entanglement is not an indi- 
cator as clear as fidelity. Indeed, considering SM modes 
of the ring, none of the entanglement derivatives seem to 
feature extrema displaying the expected scaling behaviour. 
The extrema, when present, become less pronounced with 
increasing size. The situation is more encouraging when 
the QM modes are considered. Several extrema can be sin- 
gled out that feature a scaling behaviour compatible with 
a QPT However, as it is clear from Fig. |4l a finite size 
scaling of these extrema according to Eq. (0 results in 
estimates of incompatible with each other. The fact 
that the critical point is known allows us to recognize 
that the minimum of d 2 j£ provides a satisfactory result, 
Joo = 0.262 ± 0.005 (see also Fig. E§. The rightmost 
column of Fig. [3] shows that the situation for the CCG is 
considerably simpler. The maximum in the first derivative 



of the SM entanglement provides a satisfactory estimate 
for the critical point, = 0.085 ± 0.007. All of these 
features can be understood recalling that the QPT at inte- 
ger filling belongs to different universality classes for the 
ring and CCG. This is reflected also by the fact that in all 
of the considered cases Jm approaches from different 



sides depending on the connectivity of the lattice 13 3H . In 
particular, the BKT QPT characterizing the ring has been 
recently shown to elude local measures of entanglement, 
such as the one based on SM lfl2Tl . In this respect, the 
global nature of the QM can be an explanation for the suc- 
cess of the relevant entanglement in capturing the transi- 
tion. The less elusive nature of the generic mean-field QPT 
explains the success of SM in the case of CCG. 

The last result we present is the evaluation of S (J, fi) for 
the ground state of the mean-field Hamiltonian (O. This re- 
sult is based on the (first order) analytical perturbative ex- 
pansion of the ground-state with respect to J and fi and the 
numerical evaluation of the mean-field parameter a ll23ll . 
The expected divergence of S at the superfluid side of the 
(mean-field) critical boundary is evident in the contour-plot 
in Fig. [5] An analytic perturbative expansion in a in the 
proximity of such a boundary reveals that the exponent of 
the divergence is —1. Within the Mott lobe a = 0, so 
that ground-state is independent of J and fi and 5 = 0. 
Note that the mean-field approximation "privileges" a spa- 
tial representation, resulting in a factorized ground-state, 
which prevents the identification of the QPT in terms of 
entanglement. 

According to the above results, fidelity seems to be a 
very effective tool for investigating the presence of a QPT 
and pinning down the critical point thereof. Other than the 
intuitive definition and straightforward evaluation, the key 
asset of the fidelity approach is a clear telltale of the tran- 
sition, namely a divergent behaviour of the maximum of S 
in the thermodynamic limit. Remarkably, this applies also 
for the elusive BKT transition of the ring. A further advan- 
tage of fidelity is its applicability to approximated versions 
of the model under examination, where entanglement mea- 
sures are not viable. Conversely, the more complex defini- 
tion of the considered entanglement measure allows for a 
deeper insight in the nature of the transition and underly- 
ing quantum correlations [6]. A further remarkable result 
of this work is the precise estimate of for the ID lattice. 
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FIG. 1: (color online) Mott lobe at filling v = 1. Different colors 
correspond to different sizes as specified by the color code. The 
dot signals the critical point at filling v = 1. Note that finite-size 
effects prevent the collapse of the critical boundaries. 




FIG. 2: (color online) Extrema of T and S. The vertical stripe 
and line signal Joo as reported in Refs. J20l l27ll28ll . 
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FIG. 4: (color online). Finite-size scaling of the extrema (dots) in 
the considered quantities. Solid and dashed lines denote the fits 
and the relevant estimates of Joo, as specified by the color code. 
The entanglement in the ring (CCG) is evaluated with respect to 
QM (SM). The known values of Joo are included as well. 



FIG. 3: (color online). Extrema of £ and its derivatives. The 
vertical stripes and lines signal Joo as in Fig. [2] 
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FIG. 5: Contour plot of J) for the ground-state of the mean- 
field Hamiltonian 0. 
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